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PATENT APPLICATION 
ATTORNEY DOCKET NO: 14791-501 (ARIAD-1) 

MICROARRAY-BASED SUBTRACTIVE HYBRIDIZATION 

Field of the Invention 

The present invention relates to methods for the identification and isolation of nucleic 
acid fragments. More specifically, the invention covers methods for the identification and 
isolation of non-redundant mRNAs and novel genomic sequences. 

Background of the Invention 

The elucidation of the mechanisms that dictate the normal functioning of living cells 
requires a detailed understanding of the information encoded in all of the complete genome. 
Messenger RNA (mRNA) sequences are typically used to map and sequence the genes contained 
in the genomes of different organisms. The sequence information is used to evaluate the genetic 
makeup of a particular cell or organism of interest. However, mRNAs are produced at different 
levels within different cell types and during different points in development. The distribution of 
mRNA types, their developmental and cell-type specific regulated expression, and their 
translation into protein produce the unique character of a particular cell type. 

There are currently several world-wide research efforts aimed at cloning, mapping, and 
sequencing the genomes of various organisms, including Homo sapiens. Information from these 
projects will assist in providing an understanding of how genomes result in the organisms they 
specify. Furthermore, an understanding of the molecular makeup of normal functioning cells is 
essential to the understanding of various cellular processes — including the diagnosis and 
treatment of diseases in which regulation and expression of one or more of the genes has 
changed. 

Integral to this goal is the production of libraries of cloned nucleic acids. Two different 
types of DNA libraries are typically used in the art. The first type, genomic libraries, are 
constructed by placing randomly cleaved DNA fragments of an entire genome into a suitable 



cloning vector. Although some of the clones in a genomic library contain genes or portions of 
genes, most clones contain non-coding DNA. 

The second type of DNA libraries are complementary DNA (cDNA) libraries. These 
libraries are constructed from DNA that is reverse transcribed from mRNA isolated from a 
5 source of interest; cDNA libraries primarily contain DNA that codes for genes. However, 

different species of mRNA are not equally represented in a given cell. These mRNA molecules 
are distributed into three frequency classes: (1) superprevalent (consisting of approximately 
10-15 mRNAs which, together, represent 10-20% of the total mRNA mass); (2) intermediate 
(consisting of approximately 1-2,000 mRNAs which, together, represent 40-45% of the total 
10 mRNA mass), and (3) complex (consisting of approximately 15-20,000 mRNAs which, together, 
represent 40-45% of the total mRNA mass). Davidson and Britten, Science 204\ 1052-1059 
p (1979). Differential levels of mRNAs within a cell is a significant obstacle to the identification 
?1 and sequencing of low-abundance mRNA species. In the creation of nucleic acid libraries 
'^2 suitable for sequencing, superprevalent mRNAs impede the isolation and analysis of lower 
[ft 5 abimdance mRNAs. Since the majority of clones isolated from a cDNA library will be from 
f g superprevalent and intermediate prevalent mRNAs, significant time and effort is spent 
^ resequencing previously known prevalent mRNA species, and large numbers of mRNA species 

in must be sequenced in order to isolate and sequence low-abundance mRNA species. Thus, the 

M 

^ Q rate of gene discovery from libraries is limited by the redundant nature of mRNAs present in a 
"'^0 given cell. The presence of highly abundant mRNAs also hinders the comparison of differences 
in active genes observed in different cells of related tissue types, cells in varying stages of 
development, the effect of stimuli, and differential gene expression between normally 
fimctioning and abnormal cells {e.g., cells from normal tissue compared to tumor tissues). 

One method for reducing the variation in the abundance of the individual nucleic acid 
25 molecules in a library is to produce a normaUzed library. Two approaches for generation of 
normalized libraries have been proposed. Weissman, MOL. BiOL. Med. 4: 133-143 (1987). 
These techniques include (1) hybridization to genomic DNA, in which the frequency of each 
hybridized cDNA in the normalized library is proportional to that of each corresponding gene in 
the genomic DNA, and (2) a kinetic approach that relies on the difference in annealing kinetics 
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between abundant and rare species (Galau et al, ARCH. BiOCHEM. BiOPHYS. 1 79: 584-599 
(1977). Several investigators employ the kinetic approach. For example, Soares et al. use 
single-stranded circles in their approach (see, e.g.. United States Patents 5,846,721 and 
5,830,662), while Li et al. use haptenylated nucleic acid molecules (PCT application 
5 WO 99/1 5702). An alternative approach uses reassociation of short double-stranded cDNAs. 
Ko, NUCLEIC Acids RES. 18: 5705-5711 (1990). 

Although normalization increases the chance of sequencing low-abundance nucleic acids, 
at best, the relative concentration of all mRNA species of a normalized library are within one to 
two orders of magnitude. Accordingly, the super-and intermediate-abundance nucleic acids are 
10 well represented in the library. Any attempt to randomly select and sequence clones from a 
normahzed library will result in the selection of a high percentage of previously-characterized 
f3 high abimdance nucleic acid species. 

y Therefore, a need remains in the art for a method of rapidly and efficiently identifying 

% and discarding previously-identified clones, thereby eliminating the redundancy in a population 
^35 of nucleic acid molecules. Such a method would avoid the need to continuously re-sequence 
CO previously-characterized nucleic acid fragments and would permit the rapid and efficient 
n identification and sequencing of novel genes. 

Wi 

■ 0 Summary of the Invention 

^0 The present invention relates to a highly efficient, high-throughput method for the 

20 identification and elimination of redundancy in a population of nucleic acid molecules using 
microarrays. The method comprises providing a random sample of nucleic acid fragments, 
immobiUzing the random sample of nucleic acid fragments on a microarray, hybridizing one or 
more labeled probes corresponding to previously arrayed or sequenced fragments, detecting 
fragments hybridized to the labeled probes and identifying at least one fragment not hybridized 
25 or weakly hybridized to the labeled probes; and sequencing an identified fragment that was not 
hybridized or was weakly hybridized to the labeled probes. The nucleic acid fragments may be 
RNA or DNA, and may be cloned into a vector or not. In some embodiments, the nucleic acid 
fragments are members of a cDNA or genomic library, which may be normalized or non- 
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normalized. In other embodiments, the nucleic acid fragments are PGR fragments. In many 
embodiments, the nucleic acid fragments are amplified, e.g., by PGR. 

The nucleic acid fragments are then immobilized to a solid surface, in a microarray. The 
solid surface is preferably glass. Labeled probes that correspond to previously arrayed or 
5 sequenced fragments (i.e., the subtraction probe) are hybridized to the immobilized nucleic acid 
fragments. Nucleic acid labels may be fluorescent, limiinescent, or radioactive labels, 
biotinylated, haptenated, or other chemical tags which allow for easy detection of labeled probes. 
Generally, the unhybridized probes are removed. Nucleic acid fragments that are not hybridized 
or are weakly hybridized to a labeled probe are isolated and are then pooled with the previous 
10 set of probes to generate a new, larger probe set. Usually, the newly isolated fragments are 
sequenced and their sequences compared to those found in a sequence database. 

The methods of the present invention involve a subtraction protocol that identifies and 
iJ isolates non-redundant nucleic acid fragments from a population of nucleic acid molecules. In 
> most embodiments, the protocol is reiterated, in order to create a set of fragments that becomes 
more biased toward previously uncharacterized genes with each successive round. Accordingly, 
CO with each roimd of subtraction, probes corresponding to the newly isolated fragments are labeled 
□ and added to the previous subtraction probe, and this new subtraction probe is hybridized to the 
next microarray containing randomly picked nucleic acid fragments. This procedure is repeated 
-.y several times, always adding the newly identified sequences to the previous subtraction probe. 
vSO Thus, the method allows the identification and isolation of non-redundant or minimally 

overlapping nucleic acid fragments from sources of interest and enhances the rate of novel gene 
discovery. In a preferred embodiment, the non-redundant clones that are isolated using the 
methods of the invention are identified, selected, and immobilized to a new microarray to 
produce a unigene set. 

25 Numerous applications can be envisioned for this invention. Specifically, any application 

in which the practitioner desires to enrich for sequences of interest or remove undesired nucleic 
acid fragments is amenable to the methods of the invention. A non-limiting set of uses includes: 

I. A microarray-based method for enhancing the rate of discovery of expressed 

mRNA/cDNA sequences and facilitating construction a "UniGene" set. This method 
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allows for increase novel gene discovery of expressed cDNAs and expedited construction 
of a UniGene set of expressed cDNAs. 

II. A microarray-based method for enhancing the rate of discovery of genomic sequences 
and facilitating isolation of a DNA fragments corresponding to a whole genome or 

5 subregions of interest. In this application, the method provides for increased discovery of 

genomic clones, expedited construction of a set of non redundant or minimally tiled 
genomic clones, increased discovery of clones mapping to a region of interest, expedited 
construction of a set of genomic clones in a region of interest in the genome, and 
expedited filling of gaps in genomic maps (to facilitate disease gene mapping and disease 
1 0 gene identification). 

III. A microarray-based method for enrichment and/or isolation of DNA sequences 

P (mRNA/cDNA, genomic, extrachromosomal, plasmid and all other) that are imique to a 

W population compared to another population. The invention also allows for identification 

of sequences (expressed cDNA or genomic) unique or novel to one organism versus 
-J|5 another, including nucleic acid molecules from different strains (i.e. pathogenic vs. non- 

CO pathogenic) and different species. 

IV. A microarray-based method for increasing discovery of related (or conserved) DNA 
M sequences (mRNA/cDNA, genomic, extrachromosomal, plasmid and all other). 

J5 Conversely, the method allows for identification of related sequences among closely or 

distantly related organisms. 

V. A microarray-based method for enhancing the rate of removal of undesired sequences. In 
another embodiment, the invention provides for removal of undesired DNA sequences, 
including contaminating DNA sequences and sequences closely related to those 
previously identified. 

25 VI. A microarray-based method for identifying changes in copy number (under or over 
represented) of DNA sequences (genomic, extrachromosomal, plasmid and all other) 
between different sources of nucleic acids. 



-5- 



Brief Description of the Drawings 

FIG. 1 is a flow chart depicting an overview of the gene discovery strategy. 

FIG. 2 is a representation of a microarray used in the identification of novel genes. 

FIG. 3 is a flow chart depicting the construction of a mouse calvaria unigene chip. 

5 FIG. 4 is a hne graph depicting the removal of contaminating mitochondrial genes from a 

library. 

Detailed Description of the Invention 

For convenience, the intended meaning of certain terms and phrases used herein are 
provided below. 

;-DlO "Known genes" or "known sequences" are those whose sequences can be found in a 

%^ public or private sequence database. 

Cm "Novel genes" or "novel sequences" are those which are not found in any public or 

fQ private sequence database. 

fa "Weakly" hybridizing, as defined herein, refers to a hybridization signal that has a signal- 

[15 to-noise (S/N) ratio of less than 0.5, as identified using the data analysis system of ArrayVision 

software. S/N is a measure of the background subtracted intensity of the spot signal divided by 
kB the standard deviation of the background intensity. A value of S/N less than 0.5 typically 

indicates no significant hybridization of a probe molecule to the particular DNA sample. 

"Unigene set" is a set of genes wherein each member of the set is represented once. A 
20 unigene set from a particular cell type contains one copy of each gene expressed in that cell type. 

"Complex probe" is one that contains many different DNA or RNA molecules. As used 
herein, the complex probe contains nucleic acid sequences that are complementary to known 
sequences. 

The present invention provides a highly efficient, high-throughput method for the 
25 identification and elimination of redundancy in a population of nucleic acid molecules using 
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reiterative selection involving subtractive hybridization with a complex probe pool and 
microarray technology. Randomly-chosen nucleic acid fragments chosen from a cloning library 
are microarrayed and exposed, under hybridizing conditions, to a complex labeled probe 
generated from known sequences (e.g. sequences obtained from a public database such as 
Genbank, pubEST). Only clones that do not hybridize to this first probe pool are sequenced; 
clones that hybridize to this first probe pool are not sequenced. The sequences obtained in this 
step are compared to public databases of previously-identified nucleic acid molecules (e.g. 
BLAST, pubEST, Genbank NR, Incyte). These sequences are added to the initial probe pool to 
generate a second, subtractive probe pool. This subtractive probe pool is hybridized to a second 
group of randomly-chosen clones from the library. The probe hybridizes to the nucleic acids 
complimentary to the nucleic acids of the initial, publicly known sequences, as well as to the 
nucleic acids identified in the first round of hybridization. The procedure of expanding the probe 
pool after successive rounds of hybridization is repeated until essentially all randomly-chosen 
clones hybridize to a labeled probe, indicating that all of the previously uncharacterized clones 
have been sequenced. Accordingly, the methods of the present invention eliminate the need to 
re-sequence known clones - a problem that reduces the efficiency of other methods known in the 
art. 

The first step involves isolation of randomly-chosen nucleic acid fragments from a 
population of nucleic acid molecules. These fragments may be DNA or RNA, and they may be 
cloned into vectors or not. In a preferred embodiment, the nucleic acid molecules comprise 
cDNA from a nucleic acid library of unknown sequences. In one embodiment, the library is a 
genomic library. In another embodiment, the library is a cDNA library. In either embodiment, 
the library can be a normalized library, or a non-normalized library. Alternatively, the 
population of nucleic acid fragments may not be part of a library, e.g., PGR fragments. 

In some embodiments, the nucleic acids are amplified, using PGR or other amplification 
methods. The DNA is then immobilized in a microarray on a surface, preferably glass. In some 
embodiments, the arrayed DNA comprises oligonucleotides synthesized on glass. In a preferred 
embodiment, the PGR fragments are arrayed using a microarray spotter. 




These surface-immobilized DNA molecules are then allowed to contact a labeled probe 
pool containing nucleic acid molecules of known or defined sequences under conditions which 
permit hybridization. Standard labeling protocols for nucleic acids are described, e.g., in 
Sambrook et al\ Kambara et al, BIOTECHNOLOGY 5: 816-821 (1988); Smith et al, Nuc. AciDS 
5 Res. 75: 2399-2412 (1985). Nucleic acid labels may be fluorescent, luminescent, or radioactive 
labels, biotinylated, haptenated, or other chemical tags which allow for easy detection of labeled 
probes. Fluorescent labels are advantageous for the methods described herein, as these are 
routinely used with automated instrumentation for simultaneous high throughput analysis of 
multiple samples. Metzker and Gibbs have recently disclosed a family of fluorescently tagged 
10 nucleotides based on the Cy fluorophores with improved spectral characteristics. U.S. Patent 
No. 5,728,529, incorporated herein by reference. Altemative sets of fluorophores include: the 
rhodamine based fluorophores, TARAM, ROX, JOE, and FAM; the BigDye® fluorophores 

p 

(Applied Biosystems, Inc.), the dansyl group, fluorescein and substituted fluorescein derivatives, 
n acridine derivatives, coumarin derivatives, pthalocyanines, tetramethylrhodamine, Texas Red®, 
;| 5 9-(carboxyethyl)-3-hydroxy-6-oxo-6H-xanthenes, DABCYL® and BODIPY® fluorophores 
^ (Molecular Probes, Eugene, OR). 

= Subsequent to the hybridization step, the amoimt of label is detected at each position in 

m the microarray. A key feature of the technology described here is the use of DNA microarrays. 
% See Brown and Botstein, Nat. Genet. 27: 33-37 (1999). Microarrays are an ordered array of 
'-^0 double stranded or single stranded DNA molecules positioned on a support material in a spatially 
separated organization. In contrast to filter "macroarrays", which are typically large sheets of 
nitrocellulose, microarrays position the DNA more densely packed organization such that up to 
10000 DNA molecules can be fit into a region typically 1-4 square centimeters. Microarrays 
typically use coated glass as the solid support, in contrast to the nitrocellulose-based material of 
25 fiher arrays. By having an ordered array of DNA samples, the position of each sample can be 
tracked and linked to the original sample fi:'om which the DNA on the array was generated. 
Methods and apparatus for preparing a microarray have been described. See, e.g., United States 
Patents 5,445,934 and 5,800,992, both incorporated herein by reference. 
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The process of creating a microarray involves a series of steps in which a set of bacterial 
colonies containing a cDNA library are processed to obtain an amplified fragment of purified 
DNA derived from the insert of the vector held in any one bacterial colony. The series of steps 
can be carried out either manually or through the use of robotic workstations, such as PCR 
5 machines, and liquid handling robots. In the preferred method, the steps are carried out using 
multiwell plates made up of 96 wells, 384 wells, or 1536 wells. However, other configurations 
are possible. The microarray is constructed using a robotic spotter such as the instrument 
developed by Molecular Dynamics or in the laboratory of Patrick Brown, Stanford University 
(ibid). Another key feature of the microarray process is that the data analysis is performed 
10 automatically using data analysis software, such as Array Vision. The entire process from clone 
tracking, microarray construction, to results is fiiUy integrated so that bacterial colonies to be 
processed for sequencing can be readily identified. 

The current technology for microarray construction positions up to 10,000 DNA samples 
\^ per microarray. With fiiture improvements to the technology, it will be possible to position many 
£fl5 more samples than this, perhaps 100,000 per microarray. The DNA samples on the microarray 
m are hybridized with RNA or DNA probes that have been fluorescently labeled to identify whether 

the probe sample contains a molecule that is similar or identical to the DNA sample on the 
in microarray. A key feature of the system is that a complex probe can be hybridized to the 

microarray. A complex probe is one that contains many different DNA or RNA molecules. 
'So Under the appropriate conditions, probe molecules hybridize to a DNA molecule on the 

microarray. Generally, identical or near identical sequences form productive hybrids. The 
presence of DNA-probe hybrid molecules is detected by a fluorescence detection instrument. If 
the hybridization signal is weak or non-existent at a particular DNA site, then the corresponding 
DNA or RNA molecule in the probe is absent. Current microarray instruments can hybridize up 
25 to four different fluorescent probe samples at one time. With improvements to the technology, 
more probes can be hybridized at once. 

There are several important features of the microarray-based subtraction method. For 
example, the process is accomplished by iterative rounds of microarray construction, 
hybridization, and analysis. The method reduces amount of sequencing to obtain a Unigene set 
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of a given size. It does this by removing duplicate cDNAs, samples with bacterial or 
mitochondrial DNA, and samples without any cDNA insert. 

Up until recently, DNA hybridizations were performed on nitrocellulose filters. In 
contrast to microarrays where DNA is spotted directly onto the microarray, filter arrays are 
5 generated by spotting bacterial colonies on the filters, placing the filters over a agar growth 
media, and incubating the filters under conditions that promote the bacterial colonies to grow. 
The DNA within the bacterial colonies is released by lysing the colony and treating the filters to 
fix the DNA to the filter material. The process of generating a bacterial filter array can take 
typically 2-4 days. 

10 There are a number of advantages of microarray-based subtraction versus filter array 

methods. For example, filter methods generally array bacterial colonies in which the cloned 
cDNA is contained. The colonies must be grown up over several days, lysed to release DNA and 
id fix DNA onto the filter. Hybridization to filter arrays of colonies is less reliable due to bacterial 
V debris and the low amount of DNA released from the colony. A second advantage is that the 
^JJS iterations are quicker with microarrays than with fihers. This is due to the time needed to grow 
CO colonies on the filters and prepare them for the next round of hybridization. In contrast, probing 
p of a subsequent microarray can begin less than 24 hr after analysis of an array is completed. 
Another advantage of microarrays is the ability to use fluorescently labeled probes. This 
provides for a non-radioactive method for hybridization detection. In contrast, filter 
; Jo hybridization generally uses probes labeled with radioactive phosphorus or sulfur. Microarrays 
can be hybridized with multiple probes simultaneously. In contrast, filter arrays can only be 
hybridized with one probe at a time. One of the most important advantages of microarrays is 
their reproducibility and sensitivity of hybridization signals. Typically, hybridization signals are 
higher and sensitivity is greater on microarrays versus fiUer arrays. This allows for the 
25 complexity of the probe to be higher and still get a positive hybridization signal on the 

microarray. This is a key feature for the microarray-based subtraction method. In addition, filter 
arrays often exhibit spurious background signals that are unrelated to productive hybridization 
between the probe and DNA on the filter. 
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There are a number of advantages of unigene construction by microarray-based 
subtraction versus other methods. For example, methods using selective hybridization in 
solution, as used in cDNA normalization or subtraction procedures, typically result in the loss of 
certain low abundant clones. In contrast, using the iterative microarray-based subtraction 
5 approach, all clones in the library are capable of being identified, A low abundant cDNA clone 
will be found after multiple iterative rounds of hybridization. Another advantage of the 
microarray-based approach for preparation of unigene sets of cDNAs is that the resulting cDNA 
clones are more likely to be near fiiU length. In contrast, other methods for subtraction and 
normalization often enrich short, incomplete cDNAs. 

10 A weak signal or no signal at all at a particular position of the array generally corresponds 

to the presence of DNA that has not hybridized or has weakly hybridized to sequences 
tj represented in the probe pool. Clones with a weak signal or no signal at all are sequenced. DNA 
i y sequencing methods are well known in the art. See, e.g.^ Maxam and Gilbert, Proc Natl Acad 
^; SciUSA 74: 560-564 (1977); Sanger a/., ProcNatl ACAD SciUSA 74: 5463-5467 (1977); 
£fl5 United States Patent 5,821,058, incorporated herein by reference. Automated analysis of 
f g fragments in gels or capillaries has significantly reduced the labor involved in collecting and 

processing sequence information. See, e.g., Prober et al. Science 238: 336-341 (1987); Smith 
in et al, Nature 321: 674-679 (1986); Luckey et al. Nucleic Acids Res 18: 4417-4421(1990); 
■ Q Dovichi, Electrophoresis 18: 2393-2399 (1997). The nucleic acid fragments isolated and 
'^0 sequenced in this step are added into the initial probe pool. Often, these fragments are sequenced 
first to determine the presence of non-redundant nucleic acid molecules in the microarray. Their 
sequences are then compared to public databases of previously-identified nucleic acid molecules 
(e.g., BLAST, pubEST, Genbank NR, Incyte). This second subtractive probe pool is used in the 
next round of hybridization involving a second microarray generated from a different group of 
25 randomly-chosen clones from the library. By adding the newly identified nucleic acid fragments 
to the subtractive probe pool, subsequent microarrays of randomly-chosen clones from the 
library are exposed to a larger number labeled probes. Accordingly, methods of the present 
invention eliminate the re-sequencing of known clones - a problem that reduces the efficiency of 
other art-known methods. As the number of labeled probed increases with successive iterations, 
30 the number of labeled "known" clones will increase and fewer clones will be unlabeled. 
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Eventually, with sufficient iterations, all of the clones in microarrays of randomly-chosen clones 
from a given library will be labeled, as all of the clones will have been sequenced. 

In one aspect of the present invention, a cloning library enriched in low abundance 
mRNAs is generated. In an extension of this method, the low-abimdance sequences are used in 
5 the construction of a unigene chip for use in subsequent experimentation (e.g., differential 

display experiments). An overview of the strategy for the discovery of novel genes is depicted in 
FIG. 1 . In an additional embodiment of the present invention, contaminating nucleic acids are 
removed from a cloning library. 



10 for sequences of interest or remove sequences that are not of interest. Exemplary uses include, 
but are not limited to, the following: 

J I. A MICROARRAY-BASED METHOD FOR ENHANCING THE RATE OF DISCOVERY OF 



The invention may be adapted to any purpose in which the practitioner desires to enrich 



EXPRESSED MRNA/CDNA SEQUENCES AND FACILITATING CONSTRUCTION OF A 
"UNIGENE" SET 



In one aspect, the invention provides for increased novel gene discovery and expedited 
construction of a UniGene set from expressed cDNAs. This can be accomplished by: 



(a) amplifying (by PCR or nucleic acid isolation procedures) and providing a random 
sample of nucleic acid fragments; 



(b) 



immobilizing the random nucleic acids on a solid surface in a microarray format; 



20 



(c) 



hybridizing labeled probes from a DNA source to the immobilized, microarrayed 
DNA fragments; 



(d) detecting DNA fragments hybridized to a labeled probe and identifying at least 
one fragment that does not hybridize or hybridizes weakly {i.e., subtraction); 



25 



(e) 



determining the identity of the DNA fragment by DNA sequencing, hybridization 
or other analytic approaches; and 



(f) 



reiterating steps (b) or (c) through (e) with previously identified sequences in the 
probe set in order to identify additional sequences and increase the UniGene set. 
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II. A MICROARRAY-BASED METHOD FOR ENHANCING THE RATE OF DISCOVERY OF 
GENOMIC SEQUENCES AND FACILITATING ISOLATION OF A DNA FRAGMENTS 
CORRESPONDING TO A WHOLE GENOME OR SUBREGIONS OF INTEREST 

In the case of genomic sequences, it may be desired to construct a set of "minimally 
tiled" clones, z.e., clones with the smallest amount of overlapping sequence possible. Thus, the 
invention allows for increased discovery of genomic clones, expedited construction of a set of 
non redundant or minimally tiled genomic clones and increased discovery of clones mapping to a 
gene of interest. The method also allows expedited construction of set of genomic clones in a 
region of interest in the genome. The ability to easily map clones to a region of interest in the 
genome further enables the practitioner to easily fill gaps in a genomic map, facilitating disease 
gene mapping and disease gene identification. 

The method comprises: 

(a) amplifying (by PGR or nucleic acid isolation procedures) and providing a random 
sample of genomic nucleic acid fi-agments; 

(b) immobilizing the random nucleic acids on a solid surface in a microarray format; 

(c) hybridizing labeled probes (pooled or single) from a DNA source to the 
immobilized, microarrayed DNA fragments (probes can be cDNA/mRNA or 
genomic sequences); 

(d) detecting DNA fragments which hybridize to a labeled probe; 

(e) determining the identity of the DNA fragment by DNA sequencing, hybridization 
or other analytic approaches; and 

(f) reiterating steps (b) or (c) through (e) with previously identified sequences in the 
probe set in order to identify additional sequences and increase the UniGene set. 

In particular, end regions of previously identified sequences can be used as probes to 
"walk" and identify flanking clones. 



III. A MICROARRAY-BASED METHOD FOR ENRICHMENT AND/OR ISOLATION OF DNA 
SEQUENCES (MRNA/CDNA, GENOMIC, EXTRACHROMOSOMAL, PLASMID AlVD 
ALL OTHER) THAT ARE UNIQUE TO A POPULATION COMPARED TO ANOTHER 
POPULATION. 

Another use of the invention is for the identification of sequences (expressed cDNA or 
genomic) unique or novel to one organism versus another. This includes identification of 
nucleic acids unique to one strain vs. another {Le,, pathogenic vs. non-pathogenic), as well as 
comparison of unigene sets firom closely related species and more distantly related species. 

Since nucleic acid molecules from more than one source organism are used in this 
application, the specific steps of the method are somewhat different and include: 

(a) amplifying (by PGR or nucleic acid isolation procedures) and providing a random 
sample of nucleic acid fragments; 

(b) immobilizing the random nucleic acids on a solid surface in a microarray format; 

(c) hybridizing labeled probes from source #1 and from source #2 to the immobilized, 
microarrayed DNA fragments (probes can be cDNA/mRNA or genomic 
sequences); 

(d) detecting DNA fragments which hybridize to a labeled probe from source #1 but 
not Source #2 or visa versa; and 

(e) determining the identity of the DNA fragment by DNA sequencing, hybridization 
or other analytic approaches. 

IV. A MICROARRAY-BASED METHOD FOR INCREASING DISCOVERY OF 
RELATED (OR CONSERVED) DNA SEQUENCES (MRNA/CDNA, 
GENOMIC, EXTRACHROMOSOMAL, PLASMID AND ALL OTHER) 

The converse to the above method allows discovery of related sequences (rather than 
differences) in various species, including different bacterial strains and distantly related 
organisms. Again, the method comprises: 

(a) ampHfying (by PGR or nucleic acid isolation procedures) and providing a random 
sample of nucleic acid fragments; 



(b) immobilizing the random nucleic acids on a solid surface in a microarray format; 

(c) hybridizing labeled probes (single or pooled) to the immobiUzed, microarrayed 
DNA fragments, particularly at decreased hybridization stringencies (probes can 
be cDNA/mRNA or genomic sequences); 

5 (d) detecting DNA fragments which hybridize to a labeled probe (often with a weaker 

signal); 

(e) determining the identity of the DNA fragment by DNA sequencing, hybridization 
or other analytic approaches; and 

(f) comparing DNA sequences obtained to other available DNA sequences to detect 
1 0 sequences which show homology but are not identical to other known sequences. 

° V. A MICROARRAY-BASED METHOD FOR ENHANCING THE RATE OF REMOVAL OF 
1 UNDESIRED SEQUENCES 

h The invention also provides for removal of xmdesired DNA sequences (cDNA or 

^ genomic) from any population of arrayed DNA. This may include contaminating DNA, as well 

15 as any DNA not of interest, e.g., sequences closely related to those already identified. This 

0 method comprises: 

a (a) amplifying (by PGR or nucleic acid isolation procedures) and providing a random 

5 sample of nucleic acid fragments; 

(b) immobiUzing the random nucleic acids on a solid surface in a microarray format; 

20 (c) hybridizing labeled probes, which are sequences targeted for removal, to the 

immobilized, microarrayed DNA fragments; 

(d) detecting DNA fragments hybridized to a labeled probe and identifying at least 
one fragment that does not hybridize or hybridizes weakly (Le., subtraction); 

(e) determining the identity of the DNA fragment by DNA sequencing, hybridization 
25 or other analytic approaches; and 
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(f) reiterating steps (a), (b) or (c) through (e) with previously identified sequences in 
the probe set, as deemed necessary, in order to eUminate unwanted sequences 
from the population of fragments. 



VI. A MICROARRAY-BASED METHOD FOR IDENTIFYING CHANGES IN COPY 
5 NUMBER (UNDER OR OVER REPRESENTED) OF DNA SEQUENCES (GENOMIC, 

EXTRACHROMOSOMAL, PLASMID AND ALL OTHER) BETWEEN DIFFERENT 
SOURCES OF NUCLEIC ACIDS. 

The present invention also provides a microarray-based method for identifying changes in 
copy number (under or over represented) of DNA sequences (genomic, extrachromosomal, 
10 plasmid and all other) between different sources of nucleic acids. The method comprises: 

(a) amplifying (by PGR or nucleic acid isolation procedures) and providing a random 
3 sample of nucleic acid fragments from a given source; 

^ (b) immobilizing the random nucleic acids on a solid surface in a microarray format; 

[J (c) hybridizing labeled probes (single or pooled; any type of nucleic acid - 

015 mRNA/cDNA, genomic, extrachromosomal, plasmid and all other - generally 

corresponding to the type of nucleic acid arrayed) from another source to the 

immobilized, microarrayed DNA fragments; 

sis 

'i (d) detecting DNA fragments which show absent, significantly lesser or significantly 

3 greater hybridization to a labeled probe. Such alterations in signal intensity will 

20 reflect changes in abundance; and 

(e) determining the identity of the DNA fragment(s) by DNA sequencing, 
hybridization or other analytic approaches. 

Example 1 cDNA Microarray Construction 

1 .1 cDNA Amplification 

25 A cDNA library was generated by techniques well known in the art. cDNA inserts were 

isolated from bacterial clones using Polymerase Chain Reaction (PGR). Twenty five )itl of PCR 
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mix was added to each well of a 384 well microtiter plate from a PGR master mix containing 
1000 ^il lOX PGR Buffer, 800 ^il of 2,5 mM dNTPs, 400 ^il T3 primer (5 pmoles/^1), 400 |il of 
T7 primer (5 pmoles/|il), and 10 |il recombinant Taq polymerase. The PGR plate was inoculated 
from the ovemight growth plate using a 384 pin tool to transfer approximately 1 |iL Plates were 
5 sealed with MicrosealA and placed in alpha units of MJ Research tetrad PGR machines. 

Reactions were preheated at 95 °G for 4 min, followed by 35 cycles of amplification: 45 sec at 
95°G, 1 min at 55°C, and 2.5 min at 72°G. 

Amplification of the cDNA inserts was verified by agarose gel electrophoresis in which 
2 |Lil of each PGR reaction was loaded onto a 1% agarose gel containing 500 ng/|il ethidium 
10 bromide. 500 ng of molecular weight standard (1 Kb ladder, Promega) was added to each gel for 

size determination and quantitation of amplification. Samples were electrophoresed in IX TAE 
f3 (Tris Acetate, EDTA) at 150 milliamps for 30 min. 

[:j Amplified cDNA inserts were purified from unincorporated nucleotides and primers in 

384 well glass fiber filter plates as follows: 70 |il of 5M guanidinium isothiocyanate (Sigma) 
^|5 was added to each plate using a 96 pin Gyclone liquid handler; 25 ^il of PGR reaction was 
£0 transferred from the PGR microtiter plate to the glass fiber filter plate and incubated for 2 min at 
room temperature. The fiher plate was placed on a vacuum manifold and dried. The filter plate 
J- - was washed 2 times with 70 |xl of 80% isopropanol and dried on a vacuum manifold for two min. 
^.D The purified PGR product was eluted from the glass fiber plate into a 384-well collection plate 
^|0 by addition of 50 \xl of water, incubation for 5 min at room temperature, and centrifiigation at 

3000 rpm for 5 min. Purified PGR products were lyophilized down to dryness in a speed-vac on 
high for 45 min to 1 hour. Purified products were resuspended in 30 |iil of 50% DMSO/water. 

1 .2 Immobilization of the Sample 

A maximum of 4,608 different PGR fragments were spotted onto glass slides in duplicate 
25 using a Generation III (Gen III) microarray spotter from Amersham/Molecular Dynamics. 

Twelve 384-well U bottom microtiter plates were placed in the hotel of the Gen III spotter. Up 
to 36 glass slides were placed on the spotter and spotted with DNA while the spotting humidity 
of the spotting chamber was maintained at 55%, Six slides were spotted per source plate access. 
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Prior to accessing the next set of PGR templates, the 12 pen cassette was washed as follows: 
1 sec in 0.2 M KOH, 1 sec in 95% ethanol, and 2 sec in distilled water. After spotting, slides 
were air dried for 1 hour at room temperature. The PGR fragments were cross linked to the glass 
by UV crosslinking in a Stratalinker (Stratagene) at 5000 joules. 

5 Example 2 Probe Synthesis 

The cDNA inserts for sequence verified clones were amplified in 96-well PGR plates as 
previously described for 384-well plates, with the exception that the PGR reaction volume was 
50 |al. Amplification was verified by agarose gel electrophoresis as previously described. PGR 
fragments were purified using 96-well carbon fiber filter plates and a DNA binding resin as 
10 follows: 100 |il of Wizard PGR purification resin (Promega) was added to each well of a 96-well 
filter plate using a 12 channel pipettor; the 50 |iil PGR reaction was added to the resin, and 
incubated at room temperature for 1 min; the plate was then placed on a vacuum manifold for 
1 min, until dry. Each plate was briefly washed with 200 |al of 80% isopropanol before the 
solution was removed by vacuum manifold, and the filter plate was well dried by vacuum. 
H 5 Purified PGR product was eluted by addition of 50 |il of sterile, distilled water followed by 
O centrifiigation at 3000 rpm for 5 min into a 96-well microtiter plate. The purified PGR reactions 

fcr • 

1^ were pooled and 100 |Lig of the pooled PGR reaction was purified over Qiagen PGR purification 

''S columns. The repurified PGR pool was verified by gel electrophoresis as previously described. 

^0 The concentration of the pooled PGR was determined spectrophotometrically and the 

20 concentration adjusted to 100 ng/|al using sterile, distilled water. 

The purified PGR fragment pool was used as a template for RNA synthesis in an in vitro 
transcription reaction. Purified PGR reaction (500 ng) was added to a 20 |il reaction containing 
IX transcription buffer, 0.5 mM of each rNTP, 37 mM DTT, and 10 units of T7 or T3 RNA 
polymerase. The reaction was incubated at 37°G for 90 min. The DNA template was removed 
25 by addition of 1 unit of RQl DNase and incubation at 37*'G for 30 min. The unincorporated 
nucleotides were removed by purification over an RNA Easy Golumn (Qiagen). 
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A fluorescent probe was synthesized from the in vitro synthesized RNA template in a first 
strand labeUng reaction as follows: 100 ng of RNA template was incubated with 0.5 |ig random 
hexamers in a final volume of 10 i^l for 10 min at 70°C. The reaction was chilled on ice for 
5 min, followed by the addition a reaction mix containing Ix Reverse transcription buffer 
5 (20 mM Tris pH 8.4, 50 mM KCl, 2.5 mM MgCy, 10 mM DTT, 100 |iM dGTP, dTTP, dATP, 
50 \iM dCTP, 50 \iM Cye3 dCTP or cye5 dCTP, and 200 units of Superscript 11 reverse 
transcriptase (Gibco BRL, #18089-01 1). The reaction was incubated at 42°C for 90 min. The 
RNA template was hydrolyzed by the addition of 1 ^il of 5 M NaOH and incubation at 37°C for 
10 min. The NaOH was neutralized by adding 10 ^1 of 2M MOPS (free acid). Unincorporated 
1 0 nucleotides and primers were removed by purification over GFX columns (Pharmacia). Probes 
were lyophilized down and resuspended in 30 \ii of hybridization buffer containing 50% 
formamide, 5XSSC, 0.2% SDS, IX Denhardt's, 100 \ig/m\ sahnon sperm DNA, and 1 ^g 
S oligo-dA(80). 

CJ Examples Identification of Novel Genes 

The procedure for identifying novel genes is depicted in FIG. 2. 

In 3.1 Hybridization 

sS Probes were denatured at 100°C for 10 min and added to the microarray sUde. SUdes 

were covered with glass cover slips (Coming) and incubated for 18 to 24 hours at 42°C in a 
humidified chamber. Hybridized slides were washed twice in 2X SSC containing 0.1% SDS for 
20 10 min at room temperature, followed by 2 washes in O.IX SSC containing 0.1% SDS at 42°C 
for 10 min. SUdes were dipped several times in distilled water and dried under filtered high 
pressure air. Hybridization of fluorescent probes was detected by scanning the slides in a 
Generation III (Gen III) confocal scanner (Amersham/Molecular Dynamics). 

3.2 Data Acquisition 

25 SUde images were analyzed with ArrayVision image analysis software (Imaging 

Research) for spot finding analysis, localized background determination, distribution of signal 
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intensities in a spot, and signal to noise ratios. Data was exported as a tab-delimited file and 
exported to an Oracle database. Normalization and statistical quality assessment was performed 
with a Web-based set of data analysis tools and data was analyzed with Web browser tools 
developed in the Hoechst-Ariad Genomics Center. Data was sorted based on signal to noise 
5 ratio, and DNAs which had a signal above background but below a threshold signal to noise ratio 
were further characterized by EST sequencing. 

3.3 Sequencing and Comparison to Known Libraries 

Clones which did not hybridize in step 3.1 were targeted for EST (expressed sequence 
tag) sequencing, using techniques well-known in the art. Adams a/., Science 252: 1651-1656 
10 (1991). Sequences were searched against the pubUc dbEST, mouse EST, and Lifeseq (Incyte) 

databases. All clones which were newly characterized or identified as previously-known 
i sequences were added to the probe (subtraction) pool to prevent multiple identification of the 
J same gene during subsequent hybridizations. 



fh 3.4 Reiterative Subtractive Hybridization 

f|5 The next microarray containing randomly picked clones was hybridized to a new probe 

which included all the previous clones, and in addition contained any new clones identified from 
EST sequencing. The clones identified from this methodology were used for building unigene 
1 [Ermolaeva et al, Nat. Genet. 20\ 19-23 (1998)] sets containing clones previously identified 
% from other sources (e.g., pubUc databases or Incyte) as well as novel clones which had not been 
20 previously described. 



Example 4 Selection of Low Abundance mRNAs 

A reiterative-based subtraction protocol was used where microarrays were constructed 
which contained 1500 different cDNAs randomly picked from a cDNA library. The first 
microarray was hybridized with a subtraction probe containing cDNAs encoding 64 
25 housekeeping and ribosomal genes. Non-hybridizing clones were analyzed by EST sequencing. 
Any clones which were not present in the previous subtraction probe were added to the existing 
subtraction probe and hybridized to the next microarray containing randomly picked cDNA 
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clones. This procedure was repeated 17 times, so that a total of 26,1 12 cDNAs were 
microarrayed. Based on the hybridization data, only 7,700 clones were chosen for EST 
sequencing. After clustering, 4,400 different cDNA clones were identified. This group of clones 
was highly enriched for low abundance transcripts. Additionally, by removing the redundant 
5 high abundance messenger RNAs, the microarray-based subtraction decreased the EST 
sequencing effort by 70%. 



Example 5 Construction of a Unigene Chip, use in Differential Display 

A mouse calvaria unigene chip was constructed according to the procedure outlined in 
FIG. 3. Arrays were generated of 27,648 clones from normalized mouse calvaria libraries. All 
10 27,648 clones were PGR ampUfied and arranged into 18 microarrays. As a result of subtractive 
:i hybridization (as outlined above), 7,790 cDNAs were sequenced following 18 rounds of 
p hybridization. Upon comparison of these sequences to databases containing known sequences, 
=5 4,608 clones were chosen for ampUfication and placement on a microarray. This unigene chip 

was then used in a differential gene expression experiment to identify genes modulated by BMP2 
"1 5 (bone morphogenic protein 2). 

E 

in 

Jq Example 6 Removal of Library Contamination 

J R The subtractive hybridization protocol as described by the current invention was 

demonstrated in the removal of contamination of highly abundant genes from a cloning library. 
As shown in FIG. 4, the initial library was contaminated by 35% of mitochondrial genes. After 
20 nine rounds of subtractive hybridization, this amount was reduced to 2-3%, demonstrating the 
useftilness of the current invention in improving the quality of this cloning libraries. 



Equivalents 

From the foregoing detailed description of the specific embodiments of the invention, it 
should be apparent that a unique method of gene isolation has been described. Although 
25 particular embodiments have been disclosed herein in detail, this has been done by way of 
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example for purposes of illustration only, and is not intended to be limiting with respect to the 
scope of the appended claims which follow. In particular, it is contemplated by the inventor that 
various substitutions, alterations, and modifications may be made to the invention without 
departing from the spirit and scope of the invention as defined by the claims. For instance, the 
choice of the particular source of nucleic acid, the particular label, or the particular probe pool is 
beheved to be a matter of routine for a person of ordinary skill in the art with knowledge of the 
embodiments described herein. 
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